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We present results for the lattice quark propagator in both Landau and Laplacian gauges using standard and 
improved staggered quark actions. The standard Kogut-Susskind action has errors of C(a^) while the improved 
"Asqtad" action has C?(a*), 0{a^g^) errors. This improvement is seen in the quark propagator. We demonstrate 
the application of tree-level corrections to these actions and see that Landau and Laplacian gauges produce very 
similar results. 



1. Introduction 

The quark propagator is a fundamental quan- 
tity of QCD. Though gauge dependent, it mani- 
festly displays dynamical chiral symmetry break- 
ing, contains the chiral condensate and Aqcd and 
has been used to compute the running quark mass 
(see Ref. § and references therein) . Some model 
hadron calculations rely on ansatze for the quark 
propagator, yet on the lattice we have the op- 
portunity to study it in a direct, nonperturbative 
fashion. Quark propagator studies can be compli- 
cated, however, by strong lattice artefacts ||,^. 

We are required to fix a gauge and we choose 
the Landau and the Laplacian gauges Q . We use 
Wilson glue at /3 = 5.85 (a ~ 0.125 fm) on a 16^ x 
32 lattice and six quark masses from am — 0.075 
down to 0.0125 (115 to 19 MeV). Calculations 
were done on 80 configurations. 

In the (Euclidean) continuum, Lorentz invari- 
ance allows us to decompose the full propagator 
into Dirac vector and scalar pieces 

S-\p^)=Z-\p')[i^-p + M{p')]. (1) 

Asymptotic freedom means that, as — > oo, 
S~^{p^) — > i^-p + m, (the free propagator) where 
m is the bare quark mass. 

We use the "Asqtad" quark action a fat- 
link Staggered action using three-link, five-link 
and seven-link staples to minimise flavour chang- 
ing interactions along with the three-link Naik 
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term and planar five-link Lepage term. The co- 
efficients are tadpole improved and tuned to re- 
move all tree-level 0{a?) errors. This action was 
motivated by the desire to improve flavour sym- 
metry, but has also been reported to have good 
rotational properties. 

From consideration of the tree-level forms of 
our two actions, we define the new momentum 
variables = sin(p^) and = sin(p^)[l -|- 
isin^(p^)] for the KS and Asqtad actions re- 
spectively. By considering the propagator as a 
function of q^, we ensure that the lattice quark 
propagator has the correct tree-level form and 
hopefully better approximates its continuum be- 
haviour. This is the same philosphy that has 
been used in studies of the gluon propagator (see 
Ref. Q and references therein). See also footnote 
6 in Ref. §. 

2. Laplacian Gauge 

Laplacian gauge is a nonlinear gauge fixing that 
respects rotational invariance yet is free of Gribov 
ambiguities. Although it is difficult to under- 
stand perturbatively, it is equivalent to Landau 
gauge in the asymptotic region. It is also com- 
putationally cheaper then Landau gauge. There 
is, however, more than one way of obtaining such 
a gauge fixing in SU(N). The three implementa- 
tions of Laplacian gauge fixing discussed are 

1. 9^(1) gauge (QR decomposition), used in 
Ref. §. 
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2. 9^(11) gauge, where the complex 3x3 matrix 
is projected onto SU(3) by maximising its 
trace. This will be discussed in more detail 
in an upcoming work on the gluon propa- 
gator g. 

3. 3^(111) gauge (Polar decomposition), the 
original prescription described in Ref. [Q. 

The gauge transformations employed in Lapla- 
cian gauge fixing are constructed from the low- 
est eigenvectors of the covariant lattice Laplacian 
operator. The three implementations discussed 
differ in the way that the gauge transformation 
is constructed from the above eigenvectors. We 
can think of each projection method as defining 
its own Laplacian gauge. In all cases the result- 
ing gauge is unambiguous for all configurations 
except a set of measure zero. 

3. Results 

We investigated the applicaton of tree-level cor- 
rection to the quark propagator by comparing the 
Z functions using p and q. We saw that hyper- 
cubic artefacts were small when we used p and 
they vanished almost entirely when using q. It 
is less clear which momentum variable should be 
used for the mass function, because at tree-level 
it is not multiplied by the momentum, but for 
consistency we use q here as well. In the case of 
the mass function, the choice of momentum will 
actually make little difference to our results. 

In Fig. the mass function is plotted, in Lan- 
dau gauge, for both actions with quark mass 
ma — 0.05. We see that the KS action gives 
a much larger value for M(0) than the Asqtad 
action and is slower to approach asymptotic be- 
haviour. The Asqtad action also shows slightly 
better rotational symmetry. 

The Asqtad action displays clearly better rota- 
tional symmetry in the quark Z function and, cu- 
riously, improved infrared behaviour as well. The 
Asqtad action also displays a better approach to 
asymptopia, approaching one in the ultraviolet. 
Furthermore, the relative improvement increases 
as the quark mass decreases. Comparing the mass 
function for the two actions at ma — 0.0125, the 
lowest mass studied here, the low quark mass has 
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Figure 1. Mass function for quark mass ma = 
0.05 (m ~ 77 MeV), KS action (top) and Asqtad 
action (bottom) in Landau gauge. 



introduced less noise into the propagator with the 
Asqtad action than with the KS action. 

Fig. ^ (top) shows the Z function for the Asq- 
tad action in Landau and 9^(1) gauges. Data has 
been cylinder cut Q for easier comparison. They 
are in excellent agreement in the ultraviolet but 
differ significantly in the infrared. There appears 
to be some slight difference in the Z function be- 
tween 9^(1) and 9^(11) gauges. 

Fig. H (bottom) shows the mass function for 
the Asqtad action in Landau and 9^(1) gauges. 
The two mass functions agree in the ultraviolet 
and in the infrared, but the Landau gauge mass 
function sits slightly higher in the intermediate 
region. The mass functions are nearly identical 
in 9^(1) and 9^(11) gauges. We have also found 
that Landau gauge gives slightly less anisotropy 
at this lattice spacing. 
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Figure 2. Gauge dependence of the quark Z (top) 
and mass (bottom) functions for the Asqtad ac- 
tion (ma — 0.05). Points marked with open cir- 
cles are in Landau gauge and soHd triangles are 
in 9^(1) gauge. Data has been cylinder cut. 

9^(111) performs very poorly. We found that 
many of the matrices had vanishingly small de- 
terminants (compared to numerical precision), 
which destroyed the projection onto SU(3). Prob- 
lems with (III) have also been seen in the gluon 
propagator ||]. 

A fit to each of the mass functions, including 
the chiral limit, was performed using the ansatz 



M{p) 



P 
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which is a generalisation of the one used in 
Ref. a > 1 is increasingly favoured as the 
quark mass approaches zero. We show here, in 
Fig. ^, one fit for the mass function in the chiral 
limit. 
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Figure 3. Mass function extrapolated to the chiral 
limit. Errors are Jack-knife. Fit parameters are 
c = 0.030(4), A = 870(60) MeV, mo = 0.0, a = 
1.52(23), I dof = 0.49. 

As we have simulated on only one lattice, it re- 
mains to do a thorough examination of discreti- 
sation and finite volume effects. In forthcoming 
work we will also examine the gluon propagator 
in the Laplacian gauges. 
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